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20.  ABSTRACT  (Continued) 


electric  field  were  identified,  and  a  model  of  the  most  persistent  field 
(resulting  from  the  daily  variation  of  the  magnetospheric  magnetic  field)  was 
developed.  The  model  describes  the  total  electric  field  which  includes  a  compo¬ 
nent  produced  by  the  plasma  as  it  responds  to  the  field  induced  by  the  wobbling 
dipole. 

The  natural  electric  fields  in  the  magnetosphere  have  periods  of  the  order  of 
0.25  hour,  or  more,  and  are  thus  very  low  frequency  in  character.  As  a  result 
they  penetrate  the  plasmas  of  interest  here  with  very  1 ittle  loss  of  energy  (i.e., 
have  large  skin  depths).  The  estimates  of  plasma  properties  were  performed  under 
several  simplifying  assumptions;  the  range  of  values  used  are  typical  of  the 
cases  of  interest  to  this  study.  The  dominant  feature  of  the  problem  is  the 
extremely  low  frequency  range  of  the  electromagnetic  waves.  The  normal  varia¬ 
tions  of  the  other  parameters  (such  as  the  collision  frequencies)  with  respect 
to  the  values  used  in  the  present  work  do  not  alter  the  conclusions  appreciably. 

In  addition  to  developing  the  electric  field  model  itself,  the  possible  influence 
of  this  electric  field  on  the  motion  and  structure  of  an  artificially  generated 
plasma  was  investigated.  Although  preliminary  in  nature  the  results  of  these_ 
calculations  indicate  that  the  plasma  can  be  influenced  by  this  field;  the  ExB 
drift  not  only  contributes  to  the  general  motion  of  the  plasma  but  can  also  in¬ 
fluence  internal  structuring  of  the  plasma.  High  density  plasmas  in  the  iono¬ 
sphere  are  expected  to  be  influenced  only  slightly  by  the  smaller  natural  elec¬ 
tric  field  there.  As  these  plasmas  expand  outward  along  magneti£  field  lines, 
the  penetration  of  the  field  increases  and  the  decrease  in  B  (ExB/B^  increases) 
produces  a  significant  plasma  drift  velocity.  The  penetration  of  Eis  anisotro¬ 
pic  in  nature  and  depends  on  the  detailed  structure  of  the  plasma.  The  plasma  is 
'expected  to  drift  with  respect  to  the  magnetic  field  lines,  the  variation  in 
drift  velocities  throughout  this  region  contributing  to  the  internal  structuring 
of  the  plasma. 

The  conditions  in  the  plasma  which  allow  penetration  of  the  external  electric 
field  are  also  those  under  which  the  assumption  that  magnetic  field  lines  are 
essentially  equipotential  s  becomes  poor.  I  he  breakdown  of  this  condition  al  lows 
the  less  dense  plasmas,  in  the  extended  region  of  the  magnetically  contained 
plasma  plume,  to  act  more  or  less  independently  of  the  denser,  low  altitude 
regions.  Models  of  plasma  dynamics  in  the  magnetosphere  which  depend  on  the 
"field  line  coupling"  assumption  cannot  correctly  describe  plasma  motions  in  the 
regime  of  interest  here. 

Work  was  also  performed  on  the  perturbations  to  the  ambient  fields  resulting  from 
the  presence  of  the  artificial  plasma.  The  result  is  the  evaluation  of  the  mag¬ 
netic  signatures  which  could  be  observed  by  earth  based  magnetometer  stations. 
The  results  for  a  model  high  altitude  nuclear  weapon  burst  are  compared  wi  th  the 
'available  data.  They  are  positive  and  indicate  that  this  technique  can  be  useful 
in  the  evaluation  of  the  performance  of  hydrodynamic  codes  such  as  MICE  and  MELT. 
Further  extension  of  these  techniques  to  evaluate  the  effects  of  high  altitude 
nuclear  weapon  detonations  on  the  magnetosphere  may  prove  very  useful  in  studies 
of  communication,  satellite,  and  ballistic  missile  defense  system  performance 
under  nuclear  attack.  Multiple  burst  effects  is  one  area  where  this  technique 
could  be  applied.  Also  the  techniques  used  for  quantitatively  determining 
induced  electric  fields  can  be  extended  for  use  in  the  auroral  regions  and  to 
higher  latitudes  generally  and  can  be  applied  to  the  multiburst  environment. 
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Section  1 
INTRODUCTION 


This  is  the  final  report  of  a  study  performed  at  McDonnell  Douglas 
Astronautics  Company  (MDAC)  under  contract  with  the  Defense  Nuclear  Agency 
(DNA).  MDAC  performed  the  following  tasks  enumerated  in  the  contract  work 
statement: 

1.  Provide  quantitative  information  on  the  earth's  naturally  occurring 
space  environment  as  requested  by  DNA  or  members  of  the  High 
Altitude  Nuclear  Effects  Working  Group. 

2.  Quantitatively  determine  the  total  electric  field,  in  the  presence 
of  a  plasma,  induced  by  a  naturally  occurring  time  varying  mag¬ 
netic  field. 

3.  Assess  quantitatively  the  importance  of  this  electric  field  on  the 
motions  of  plasmas  artificially  injected  into  the  geomagnetic 
field. 

During  the  contract  effort  representatives  from  MDAC  also  attended 
several  phenomenology  working  group  meetings  and  met  several  times  with  the 
DNA  contract  monitor.  Detailed  exchanges  of  technical  information  were 
also  held  between  the  MDAC  group  working  on  this  contract  with  several  members 
of  the  Mission  Research  Corporation  (MRC)  organization  associated  with  the 
running  of  the  hydrodynamic  codes  MICE  and  MELT. 

From  these  discussions  it  became  apparent  that  any  technique  which 
could  be  used  to  verify  the  validity  of  these  high  altitude  hydrodynamic 
calculations  would  be  of  interest  to  DNA.  As  a  result,  work  not  included  in 
the  contract  work  statement  was  performed  in  an  attempt  to  quantitatively 
model  the  magnetic  signature  of  a  high  altitude  burst. 

Techniques  developed  at  MDAC,  used  in  generating  models  of  the 
magnetospheric  magnetic  field,  were  applied  to  the  internal  current  systems 
of  high  altitude  artificially  produced  plasmas  to  calculate  their  magnetic 
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signatures  at  any  location  in  and  around  the  plasmas.  Preliminary  results 
are  included  in  this  report.  The  plasma  plume  model  and  internal  current 
distributions  used  in  these  calculations  were  obtained  from  MRC  and  result 
from  calculations  performed  by  them  with  the  MICE/MELT  code.  We  recommend 
that  similar  calculations  be  performed  by  MRC  for  actual  weapon  test  events 
and  the  MDAC  procedures  can  then  be  used  to  compare  model  results  with 
the  observed  perturbations  of  the  geomagnetic  field.  This  is  a  straight¬ 
forward  and  available  way  to  test  the  accuracy  of  the  late  time  hydrodynamic 
codes. 

Included  in  this  report  is  a  brief  description  of  the  magneto¬ 
sphere,  emphasizing  those  features  which  may  have  a  direct  bearing  on 
artificially  produced,  high  altitude  plasmas.  A  brief  description  of  a 
time  dependent  model  of  the  earth's  external  magnetic  field  developed  at 
MDAC  is  also  provided.  The  induced  electric  field  associated  with  the 
model  magnetic  field  is  also  described.  It  is  given  in  terms  of  the  time 
varying  magnetic  vector  potential  and  is  used  to  quantitatively  assess 
the  penetration  and  persistance  of  the  naturally  occurring  electric  field 
in  a  dense  plasma.  The  techniques  used  to  generate  the  magnetic  induction 
vector,  F,  in  the  magnetic  field  model  and  the  vector  magnetic  potential, 

A,  used  to  calculate  the  electric  field  are  both  described.  These  techniques 
were  modified  and  applied  to  the  calculations  of  magnetic  signatures  of 
high  altitude  artificial  plasma. 

The  magnetospheric  electric  field  model  presented  in  this  report 
represents  only  those  fields  which  results  from  the  daily  wobble  of  the 
geomagnetic  dipole.  Other  naturally  occurring  electric  fields  are  much 
larger  in  magnitude  but  are  more  limited  in  extent  and/or  duration.  The 
particular  choice  of  induced  electric  field  made  for  this  study  was  pri¬ 
marily  dictated  by  the  availability  of  a  quantitative  model  of  both  F  and 
A"  caused  by  the  wobbling  dipole. 

The  latest  MDAC  magnetic  field  model  (also  described  here)  is 
sufficiently  flexible  to  include  fluctuations  in  the  magnetic  field 
resulting  from  any  variations  in  the  magnetospheric  current  system.  It 
is  capable  of  modeling  magnetic  substorm  activity  and  any  other  perturba¬ 
tions  of  the  magnetospheric  fields  resulting  from  changes  in  the  strengths 
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and  topologies  of  the  current  systems.  High  altitude  nuclear  weapon 
detonations  (such  as  those  occurring  at  satellite  altitudes)  can  also 
perturb  the  magnetospheric  current  systems  and  be  responsible  for  large 
scale  fluctuations  in  the  earth's  magnetic  field.  The  energy  deposited 
in  this  system  may  produce  subsequent  substorm  activity.  The  techniques, 
used  to  model  the  general  time  varying  behavior  of  the  naturally  occurring 
electromagnetic  field,  can  be  extended  for  use  in  determining  the  total 
environment  in  a  multiburst  scenario.  (It  has  already  been  used  to  study 
burst  magnetic  signatures  as  is  discussed  in  this  report.) 

In  the  sections  of  this  report  dedicated  to  assessing  the 
interaction  of  artificially  produced  plasmas  with  the  natural  magneto- 
spheric  fields,  several  aspects  of  this  interaction  are  investigated. 

Of  principal  concern  to  this  study  was  the  interaction  of  the  electric 
field  with  the  plasma  and  the  resulting  plasma  structuring  and  motion. 

This  was  investigated  for  the  case  of  the  electric  field  induced  by  the 
diurnal  variation  of  the  magnetic  field.  Various  regimes  for  this  plasma 
are  identified  where  even  this  rather  weak  electric  field  can  have  a 
significant  effect  on  plasma  structuring  and  motion.  Fields  which  result 
from  other  naturally  occurring  sources  are  also  expected  to  be  significant 
in  this  respect  but  have  not  yet  been  quantitatively  examined. 
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Section  2 

OBSERVED  ELECTROMAGNETIC  ENVIRONMENT 


In  order  to  properly  describe  the  structure  of  the  plume  and  the 
motion  of  low  energy  plasma  during  late  times,  it  is  necessary  to  understand 
the  natural  plume  environment  in  the  earth's  ionosphere  and  magnetosphere. 

In  this  section  the  particle  and  field  properties  of  these  regions  are 
briefly  described  with  attention  given  to  the  electromagnetic  field. 

The  earth's  magnetosphere  is  formed  and  maintained  by  the  inter¬ 
action  of  the  earth's  magnetic  field  with  the  solar  wind.  The  solar  wind 
is  a  neutral  plasma  whose  velocity  (in  an  earth  reference  frame)  varies 
from  about  300  to  1000  km/sec  and  whose  density  of  ion-electron  pairs 
varies  from  0.1  to  several  hundred  per  cubic  centimeter.  The  solar  wind 
carries  with  it  a  magnetic  field  which  is  seen,  when  transformed  into  an 
earth  reference  system,  as  both  a  magnetic  and  electric  field.  Because  of 
the  highly  variable  nature  of  the  solar  wind,  the  magnetosphere  exhibits 
both  "quiet"  and  "disturbed"  topologies. 

The  total  magnetospheric  magnetic  field  is  important  in  the  quan¬ 
titative  description  of  the  nuclear  burst  in  several  ways.  First  the  burst 
particles  must  expand  against  the  pressure  of  the  magnetic  field.  Later  the 
magnetic  field  topology  determines  the  range  over  which  the  low  energy 
charged  particles  move  and  also  the  location  of  the  conjugate  region. 

The  earth's  main  magnetic  field  exhibits  only  a  secular  variation 
and  for  the  nuclear  plume  problem  can  be  represented  as  constant  in  time. 
Also,  for  most  applications  it  suffices  to  represent  the  main  field  as  a 
dipole  with  a  pole  strength  of  0.314  gauss  times  R£  (RE  =  earth  radii).  This 
dipole  is  placed  at  the  center  of  the  earth  with  an  inclination  of  11.7°  to 
the  earth's  rotation  axes  and  lying  in  the  meridian  plane  69°  west  of 
Greenwich.  For  the  accurate  calculation  of  the  field  in  the  vicinity  of  the 
earth  it  is  sometimes  necessary  to  use  a  more  complex  representation  of  the 
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main  field.  Such  descriptions  are  available  as  spherical  harmonic  expansions.  , 

Such  representations  may  be  necessary  to  accurately  describe  particle  and  < 

field  behavior  at  a  particular  point  in  space  and  also  to  determine  the 
motions  of  charged  particles  and  locations  of  magnetic  conjugate  regions. 

The  formation  of  the  magnetosphere  by  the  dipole  -  solar  wind 
interaction  is  caused  mainly  by  the  magnetopause  current  system.  These 
currents  are  formed  by  the  deflections  of  the  solar  wind  plasma  by  the  geo¬ 
magnetic  dipole.  (Since  the  solar  wind  electrons  and  protons  are  deflected 
in  opposite  directions,  they  form  a  current  system.)  These  currents  grossly 
distort  the  magnetic  field  within  the  magnetosphere  and  simultaneously  limit 
the  penetration  of  the  solar  wind.  The  region  over  which  these  currents 
flow  is  called  the  magnetopause  -  the  boundary  of  the  magnetosphere.  Thus, 
the  magnetosphere  may  be  thought  of  as  that  region  of  space  where  the  earth's 
magnetic  field  persists. 

Two  other  current  systems  flow  in  the  magnetosphere  which  also 
grossly  influence  magnetic  field  topology.  In  the  inner  magnetosphere  where 
the  field  geometry  remains  essentially  dipolar,  charged  particles  can  remain 
"trapped"  for  long  periods  of  time.  Since  the  particle  motions  are  charge 
dependent,  a  torroidal  current  system,  the  so-called  ring  current,  is  formed. 

In  the  anti  solar  region  (the  "tail"  of  the  magnetosphere)  currents  flow 
across  the  equatorial  region  much  the  same  way  they  did  in  the  ring  current 
except  that  they  are  not  permanently  trapped  but  flow  in  and  out  of  the  tail 
along  its  flanks. 

In  order  to  properly  describe  the  interaction  of  the  plume  during 
late  times,  it  is  necessary  to  describe  the  total  magnetic  field  produced 
by  the  earth's  main  field  and  the  fields  from  these  magnetospheric  current 
systems.  Near  the  earth's  surface,  of  course,  the  main  field  dominates  and 
it  may  not  be  necessary  to  consider  the  magnetospheric  contribution. 

However,  at  geosynchronous  altitude,  the  magnetospheric  currents  produce  a 
large  portion  of  the  total  field.  For  L  values  larger  than  4  it  is  generally 
necessary  to  consider  the  magnetospheric  contributions.  (The  L  value  is  the 

4 
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distance  of  the  farthest  extent  of  the  field  line  from  the  center  of  the 
earth  measured  in  earth  radii.  Thus,  the  L  value  of  field  lines  through 
geosynchronous  orbit  is  about  6.6.) 

Although  the  main  field  can  be  considered  as  constant  with  time, 
the  major  magnetospheric  currents  systems  are  highly  variable  as  a  result 
of  their  control  by  the  solar  wind.  The  magnetopause  currents  for  example 
typically  cross  the  earth  sunline  approximately  10  up  stream  of  the 
earth's  center  during  quiet  conditions  (typically  associated  with  a  steady 
low  speed  solar  wind).  However,  the  magnetopause  is  on  occasion  observed 
to  be  inside  geosynchronous  orbit.  During  such  periods  the  magnetic  field 
is  considerably  distorted  and  the  behavior  of  both  naturally  occurring  and 
artificially  injected  charged  particles  is  considerably  different  than 
during  quiet  conditions. 

Both  naturally  occurring  and  artificially  produced  time  variations 
in  the  magnetic  field  cause  currents  to  flow  in  the  earth's  crust  which  tend 
to  shield  the  interior  of  the  earth  from  the  magnetic  variations.  These 
telluric  currents  also  distort  the  magnetic  field  above  the  earth's  surface. 
They  have  been  studied  in  considerable  detail  with  relation  to  naturally 
occurring  magnetic  variations.  The  distortion  of  the  local  magnetic  field 
produced  by  telluric  currents  should  be  included  in  any  study  of  the  time 
varying  magnetic  field  produced  by  a  low  altitude  nuclear  burst.  An  accurate 
computation  of  the  field  at  the  location  of  the  ground  based  magnetometers 
must  include  the  effects  of  these  telluric  currents. 

In  addition  to  the  problems  this  variability  creates  for  the 
accurate  description  of  the  magnetic  field,  the  time  varying  nature  of 
magnetic  field  produces  an  induced  electric  field  that  can  strongly  influence 
the  behavior  of  low  energy  charged  particles.  The  subject  of  induced  elec¬ 
tric  fields  is  discussed  in  the  paragraphs  below. 

Low  energy  charged  particles  are  influenced  by  electric  fields  in 
the  magnetosphere.  The  magnetospheric  electric  field  has  many  sources 
(Table  2-1  lists  some  of  the  sources  and  gives  the  approximate  electric  field 
strength).  These  include  electrostatic  fields  produced  by  charge  separation 
along  field  lines.  This  phenomenon  appears  to  occur  in  the  auroral  region 
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and  is  the  source  of  acceleration  for  the  auroral  primary  particles.  It 
also  acts  to  accelerate  heavy  ions  into  geosynchronous  orbit  from  the 
ionosphere.  An  electric  field  associated  with  the  convection  of  low  energy 
plasma  in  the  magnetosphere  that  is  produced  by  a  "viscous  interaction" 
between  the  solar  wind  and  magnetospheric  plasma  has  been  widely  discussed 
but  is  not  well  understood  quantitatively. 


Table  2-1 

ESTIMATES  OF  INDUCED  ELECTRIC  FIELDS 


B 

t 

£s 

E 

Source 

(gamma) 

(hours) 

(RE} 

(volts/m) 

Daily  wobble  of  dipole 
plus  image 

30 

24.0 

10.0 

0.0005 

Field  aligned  currents 

250 

0.25 

0.1 

0.004 

Tail  field  collapse 

50 

0.25 

10.0 

0.09 

Sudden  commencement 
magnetic  storm 

1500 

0.25 

10.0 

2.5 

Electric  fields  induced  by  time  variations  in  the  magnetospheric 
magnetic  fields,  however,  have  been  quantitatively  represented  and  are  quite 
well  understood.  One  of  the  key  goals  of  this  contract  has  been  to  under¬ 
stand  the  ways  in  which  these  naturally  occurring  induced  electric  fields 
influence  gross  motions  and  structuring  in  the  plume  of  a  nuclear  burst 
during  late  times.  The  magnetospheric  electric  field  generally  can  be 
characterized  by  its  variability.  During  "quiet  times"  the  magnitude  of  the 
fluctuations  in  the  electric  field  are  typically  as  large  as  the  field  itself 
from  the  DC  field  to  frequencies  of  several  hertz.  During  disturbed 
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conditions  the  total  electric  field  strength  can  increase  dramatically  and 
the  fluctuations  are  a  smaller  fraction  of  the  total  magnitude. 

The  various  sources  of  this  induction  electric  field  and  its 
strength  are  discussed  elsewhere  in  this  report.  In  this  contract  effort 
we  have  dealt  almost  exclusively  with  the  electric  field  induced  by  the 
daily  wobble  of  the  earth's  dipole  axes.  Although  other  naturally  occurring 
induced  fields  are  much  larger  in  magnitude,  this  particular  source  was 
studied  because  it  is  always  present  and  it  is  well  understood. 
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ELECTROMAGNETIC  FIELDS 


Models  of  the  quiet  time  magnetos pheric  magnetic  and  electric 
fields  are  described  in  this  section;  also  described  are  the  possible  exten¬ 
sions  of  these  models  to  include  effects  due  to  a  disturbed  magnetosphere. 

The  model  of  the  earth's  external  magnetic  field  was  developed 
under  other  contracts  and  are  fully  documented  elsewhere  (Olson  and  Pfitzer, 
1977).  A  brief  description  of  the  magnetospheric  magnetic  field  is  provided 
here  in  order  to  illustrate  an  example  of  naturally  occurring  time  dependence 
of  the  magnetic  field  which  has  associated  with  it  an  induced  electric  field 
which  is  described  insubsequent  sections.  (The  same  techniques  used  to 
obtain  this  magnetic  field  were  later  modified  to  allow  calculations  of 
magnetic  signatures  of  artificial  plasma  clouds  in  the  magnetosphere.)  These 
techniques  are  quite  general  and  can  be  applied  to  other  plasma-field  prob¬ 
lems  such  as  the  nuclear  weapon,  multiple  burst  considered  elsewhere  in  this 
report. 

3.1  THE  MAGNETIC  FIELD  MODEL 

In  order  to  quantitatively  determine  the  magnetic  field  in  near 
earth  orbit,  it  is  necessary  to  first  obtain  quantitative  representations  of 
the  magnetospheric  current  systems.  They  include  the  magnetopause  current 
system,  formed  by  the  interaction  of  the  solar  wind  particles  with  the  total 
magnetospheric  magnetic  filed,  the  tail  current  system  produced  by  particles 
flowing  across  the  tail  of  the  earth's  magnetosphere  and  returning  on  its 
boundary,  and  the  quiet  time  ring  current  system  produced  by  the  motion  of 
charged  particles  moving  within  the  magnetosphere.  The  magnetopause  current 
is  formed  by  particles  which  interact  with  the  earth's  magnetic  field  and 
then  return  immediately  to  the  magnetosheath.  These  particles  interact  with 
the  current  system  only  for  fractions  of  seconds  as  they  are  deflected  by 


the  magnetospheric  magnetic  field.  The  particles  which  form  the  tail 
currents  remain  part  of  that  current  system  for  as  long  as  several  days  as 
they  make  their  journey  across  the  tail.  The  quiet  time  ring  current  system 
is  formed  by  particles  trapped  along  magnetic  field  lines  in  the  inner 
magnetosphere.  These  particles  may  contribute  to  the  currents  for  as  long 
as  weeks  or  even  months.  Although  the  magnetopause  currents  flow  only  on  a 
surface  (on  the  magnetopause  itself)  both  the  tail  current  system  and  the 
ring  current  particles  are  distributed  over  a  large  region  of  space  in  the 
magnetosphere.  Thus,  the  ring  and  tail  currents  are  initially  represented 
in  the  form  of  wire  loops.  In  the  inner  magnetosphere  there  is  a  nest  of 
three  sets  of  wire  loops  used  to  represent  the  ring  current.  The  tail 
currents  flow  almost  linearly  across  the  upper  and  lower  boundaries  of  the 
plasmasheet  in  the  distant  tail  while  close  to  the  earth  they  curve  around 
such  that  they  are  in  close  proximity  to  the  ring  current  loops.  The  return 
path  of  the  tail  currents  was  constructed  to  flow  approximately  on  the 
boundary  of  the  magnetopause  with  the  shape  similar  to  that  used  in  the  deter¬ 
mination  of  the  magnetopause  currents.  A  three-dimensional  view  of  these 
current  systems  is  shown  in  Figure  3-1. 

The  procedure  for  determining  the  appropriate  values  of  these 
currents  was  as  follows.  The  general  form  of  the  currents  as  suggested  by 
the  observations  were  used  to  select  the  positions  of  "wires"  and  the  strength 
of  the  current  flowing  through  each  of  them.  Integration  over  these  currents 
then  yielded  values  of  the  field  from  the  distributed  currents  that  were 
added  to  the  magnetopause  magnetic  field  (which  had  been  previously  deter¬ 
mined)  and  the  earth's  main  field  and  compared  with  observations.  Differences 
in  the  observed  and  computed  values  of  the  field  were  noted  and  changes  made 
in  the  positions  and  strengths  of  the  currents  in  order  to  decrease  the 
differences.  The  first  check  on  the  model  is  to  compare  it  with  the  observed 
magnitude  of  the  external  field.  Other  tests  of  the  accuracy  of  the  model 
are  given  below. 

The  total  magnetospheric  magnetic  field  can  then  be  found  by 
integrating  over  these  external  currents  and  adding  the  field  produced  by 
the  earth's  main  field.  However,  this  is  a  very  time  consuming  process.  For 
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the  model  to  be  useful  it  must  be  capable  of  returning  the  total  magnetospheric 
magnetic  field  at  any  position  in  the  magnetosphere  with  a  minimal  expenditure 
of  computer  time.  Therefore,  the  values  of  the  magnetic  field  were  found 
at  many  points  by  direct  integration  over  the  distributed  and  magnetopause 
currents.  These  values  of  the  field  were  then  input  to  a  least  squares  best 
fit  computer  program  to  determine  coefficients  of  the  polynominal  series  to  be 
used  in  the  analytic  field  representation.  Previous  representations  of  the 
magnetopause  magnetic  field  have  been  given  in  the  terms  of  the  scalar 
potential.  Since  the  field  from  the  distributed  currents  must  be  fit  in  the 
region  of  the  source  currents,  a  scalar  potential  can  no  longer  be  used  since 
VxB  must  be  permitted  to  be  nonzero.  Another  constraint  on  the  form  of  the 
series  representation  is  dictated  by  the  rather  sharp  structuring  that  the 
field  from  the  distributed  currents  exhibits  in  the  inner  magnetosphere.  Thus, 
the  series  must  simultaneously  fit  a  very  structured  region  near  the  earth  and 
the  weak  field  region  in  the  distant  tail.  A  standard  one-dimensional  ortho¬ 
normal  least  squares  fitting  program  was  generalized  for  use  to  find  coeffi¬ 
cients  in  n  variables.  In  order  to  fit  both  the  structure  of  the  inner 
magnetosphere  and  the  tail  field  topology,  it  was  necessary  to  include  terms 
of  up  to  sixth  order  in  a  power  series  plus  similar  terms  multiplied  by  an 
exponential.  A  final  requirement  for  the  model  is  that  the  field  must  vary 
smoothly  between  the  points  where  it  was  computed  for  input  to  the  fitting 
program.  Calculations  of  both  the  distributed  magnetic  field  and  the  magneto¬ 
pause  magnetic  field  show  that  the  model  fields  do  vary  smoothly  over  the 
region  where  they  are  defined.  This  finding  is  particularly  important  for 
the  distributed  magnetic  field  since  "wires"  were  used  initially  to  represent 
source  currents.  The  form  of  the  series  representation  for  the  distributed 
magnetic  field  and  the  coefficients  have  been  published  by  Olson  and  Pfitzer 
(1974). 

Field  lines  representing  the  total  quiet  time  magnetospheric  magnetic 
field  are  shown  in  the  noon-midnight  meridian  plane  for  two  angles  of  incidence 
of  the  solar  wind  on  the  dipole  axis  in  Figure  3-2.  It  is  this  daily  variation 
in  tilt  angle  and  the  associated  3B/3t  that  induces  an  electric  field  throughout 
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3-2. 


the  magnetosphere.  The  quantitative  modeling  of  this  electric  field  and 
induced  electric  fields  in  general  is  discussed  in  the  next  section.  Note 
that  the  field  line  from  76°  magnetic  latitude  cross  the  equator  at  midnight 
farther  from  the  earth  than  the  82°  lines  in  previous  models.  It  is  this 
feature,  the  extension  of  the  near  earth  field  lines  produced  by  the  inclusion 
of  the  distributed  currents  that  makes  this  model  capable  of  reproducing  many 
particle  and  field  observations  that  previously  were  not  properly  modeled. 

The  accuracy  of  the  model  has  been  tested  by  comparing  model 
calculations  with  particle  and  field  observations.  Optical  tracking  data 
(Adamson,  et  al . ,  1973)  from  the  NASA-MPE  Barium  Cloud  Experiment  has  shown 
that  the  field  lines  in  the  inner  magnetosphere  are  more  elongated  than  the 
lines  calculated  from  models  consisting  of  main  field  and  earlier  external 
field  representations  (Barrish  and  Roederer,  1973).  A  computer  program  that 
combines  any  series  expansion  of  the  main  field  with  various  representations 
of  the  magnetospheric  field  used  the  IGS75  main  field  expansion  and  the 
present  model  to  generate  field  lines  and  compare  them  with  the  optical  track 
of  the  barium  clouds  (see  Figure  3-3).  The  model  field  fits  the  observations 
quite  accurately  because  it  includes  the  depressed  inner  magnetosphere  field 
feature  produced  by  the  quiet  time  ring  current.  The  earth  intercept  of  field 
lines  from  geosynchronous  orbit  have  been  determined  for  various  local  times 
(see  Figure  3-4).  The  foot  of  the  field  line  lies  between  65°  and  66° 
magnetic  latitude  depending  on  the  time  of  day.  This  is  almost  3°  lower 
than  the  values  given  by  previous  boundary  models. 

3.2  THE  ELECTRIC  FIELD  MODEL 

Current  understanding  of  the  magnetospheric  electric  field  is 
incomplete.  Although  data  on  this  electric  field  are  becoming  increasing 
available,  there  is  currently  no  quantitative  description  of  the  various 
sources  of  the  magnetospheric  electric  field.  The  one  source  of  these  fields 
which  is  relatively  well  understood  is  the  time  variation  in  the  magnetic 
field.  A  model  of  the  quiet  time  tilt  dependent  magnetic  field  was  examined 
in  the  preceding  section,  and  its  extension  to  include  distributed  magneto¬ 
spheric  effects  is  discussed  in  the  following  section  (Section  3.3).  These 
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•  OPTICAL  TRACKING  DATA 
(ADAMSON,  ET.  AL,  1973) 
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Figure  3-4.  Earth  Intercept  of  Field  Lines  for  Geosynchronous  Orbit. 


models  consider  the  time  variations  in  all  the  major  magnetospheric  current 
systems  as  sources  of  the  total  external  magnetic  field. 

The  electric  field  model  described  in  this  section  results  from  the 
diurnal  fluctuations  of  the  magnetic  field  model  caused  by  the  solar  wind 
plasma  -  geomagnetic  dipole  field  interaction.  It  is  emphasized  that  the 
procedure  used  to  obtain  this  electric  field  can  equally  well  include  any 
electric  fields  induced  by  temporal  fluctuations  in  the  magnetic  field. 

The  procedure  described  in  Section  3.1  described  models  of  the 
magnetospheric  current  systems;  the  Biot-Savart  law  is  applied  to  these 
currents  to  yield  the  magnetic  field,  thus 


B(F)  .  (  3(P)  x  _ir=rO-  dV'. 

4tt  J  |r-r'| 

Because  of  the  complexity  of  this  current  system,  the  integral  is  performed 
numerically.  Similarly,  and  with  very  little  additional  cost  in  computer 
time,  the  vector  potential  for  this  field  can  also  be  obtained. 
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Since  the  magnetospheric  current  system  model  includes  a  dependence  on  the 
daily  precession  of  the  earth's  dipole  moment,  this  particular  time  depen¬ 
dence  of  the  magnetic  field  and  the  associated  vector  potential  are  properly 
modeled.  The  induced  component  of  the  magnetospheric  electric  field  resulting 
from  this  dependence  of  the  magnetic  field  can  be  obtained  from  the  vector 
potential 
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Because  of  the  presence  of  charged  particles  (low  density  plasma)  throughout 
the  magnetosphere  (this  medium  is  the  source  of  the  current  system),  local 
electric  fields  resulting  from  this  medium  must  also  be  considered. 

Indeed,  the  total  electric  field  at  a  given  location  will  not 
only  include  the  contribution  from  the  local  time  dependence  of  the  magnetic 
field,  but  also  contributions  from  other  regions  of  the  magnetosphere 
transmitted  to  this  location  by  the  medium,  to  the  extent  that  it  is  conductive. 
Thus,  the  total  electric  field, 

ET  •  ■  -  3t  ' 

at  a  given  point  in  the  magnetosphere  includes  contributions  from  the 
induced  field  from  other  regions  of  the  magnetosphere,  linked  to  the  given 
location  by  the  conductivity  of  the  medium,  as  well  as  any  contributions 
from  location  by  the  conductivity  of  the  medium,  as  well  as  any  contributions 
from  local  charge  separation  in  the  medium.  This  dependence  is  contained  in 
the  scalar  potential  component,  E<.  =  of  the  total  field. 

A  procedure  developed  by  Hones  and  Bergeson  (1965)  was  applied  to 
evaluate  the  scalar  potential  component  of  the  electric  field.  It  is  assumed 
that  in  the  region  of  interest,  charged  particle  densities  are  high  enough 
to  assure  sufficient  conductivities  along  the  magnetic  field  lines  such  that 
Ft-B  =  0. 

Then 

d£  _  9An_ 

ds  3t  ’ 

where  d<{>/ds  is  the  scalar  potential  gradient  along  the  field  line,  and  9A„/3t 
is  the  component  of  fTj  parallel  to  the  magnetic  field  line.  This  equation 
can  be  integrated  along  the  field  line  to  give  the  scalar  potential  at  a 
given  point  in  the  magnetosphere.  Thus, 
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where  the  line  integral  is  along  the  field  line,  and  <j>g  is  the  potential 
at  some  appropriate  boundary. 

In  the  modeling  of  a  realistic  total  electric  field,  care  must 
be  taken  in  the  selection  of  the  proper  boundary  condition  as  well  as  the 
value  of  3A„/3t  to  be  used  in  the  integral.  The  presence  of  the  medium 
complicates  matters  and  careful  estimates  of  the  charged  particle  densities 
along  the  integration  path  must  be  made.  If  conductivities  along  the  field 
line  are  sufficiently  high,  the  induced  electric  field  will  propagate  through 
very  short  distances  and  only  local  contributions  of  Ej  will  be  significant. 
At  lower  densities,  the  integral  may  be  extended  along  the  magnetic  field 
lines  but  then  the  question  of  the  appropriate  boundary  condition  must  be 
considered.  The  various  options  available  for  the  selection  of  3A„/3t  are 
summarized  in  Table  3-1. 


Table  3-1 
VALUE  OF  %- 


3^total  _ 
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if  E*B  =  0  ,  then  =  -rr-  (Hones  and 
ds  at  „ 

Bergeson; 

Ionized  Media 

Assumed) 

by  energy  dissipation  A.  dcp/ds  =  3A„/3t| local 

B.  d<t>/ds  'v  0  due  to  high 
conductivity 

if  E»B  f  0,  then  d4>/ds  =  ? 
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The  simple  case  of  time  independent  B  is  included  here  only  for 
completeness,  since  it  does  not  represent  a  very  interesting  case.  It  is 
consistent,  however,  with  the  assumption  of  magnetic  field  lines  being 
equipotentials  in  an  ionized  medium.  The  second  case  represents  the  Hones 
and  Bergeson  procedure  (discussed  above)  which  is  appropriate  in  regions  where 
the  conductivity  along  field  lines  is  high  enough  so  that  E*B  =  0.  The 
principal  questions  in  applying  this  procedure  is  the  selection  of  an 
appropriate  boundary  condition  and  the  sources  of  magnetic  field  which  con¬ 
tribute  to  9A„/3t.  In  the  limit  of  very  high  conductivities,  3A„/9t  approaches 
zero  through  energy  dissipation,  and  only  local  sources  of  the  magnetic  field 
contribute.  Case  3  illustrates  this  condition;  again  the  magnetic  field  lines 
will  either  approximate  equipotentials  or  only  local  sources  of  3A„/3t  will 
contribute  to  E  parallel  to  the  magnetic  field  line. 

Case  4,  in  Table  3-1,  represents  the  conditions  present  in  areas 
where  the  conductivity  along  field  lines  is  so  low  that  F-B  f  0.  The  pro¬ 
cedures  discussed  in  this  reportare  no  longer  appropriate  for  this  case  and  the 
only  available  approach  is  the  solution  of  the  wave  equation.  This  could  be 
done  if  a  particularly  convenient  boundary  condition  for  the  regions  of 
interest  could  be  postulated. 

Options  for  the  selection  of  the  appropriate  boundary  condition  for 
<p  include  the  following: 

1.  Experimentally  determined. 

2.  Constant  over  an  appropriate  surface. 

3.  Uniformally  magnetized  rotating  sphere  (Hones  and  Bergeson, 

1965), 

V  =  U-A  , 

where  V  =  the  potential  at  point  r, 

U  =  w  x  7, 

A  =  y  x  r/R^ 
w  =  angular  velocity 
y  =  magnetic  moment 
Rs=  radius  of  sphere. 

4.  A  self-consistent  rigorous  determination  of  <p  and  A. 
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In  the  present  work,  options  2  and  3  were  employed  to  evaluate 
the  sensitivity  of  the  electric  field  model  to  both  the  choice  of  the 
boundary  condition  and  the  nonlocal  contributions  to  the  induced  electric 
field.  Some  results  of  this  investigation  are  illustrated  in  Figures  3-5 
and  3-6.  In  both  figures,  equipotential  contours  are  plotted  in  the  earth's 
equatorial  plane.  The  positive  x-direction  in  both  plots  is  through  the 
Greenwish  meridian.  The  calculations  were  done  in  an  inertial  coordinate 
system  in  which  the  earth  is  rotating,  the  induced  electric  field  (i.e., 

3A„/3t)  is  due  to  both  the  dipole  field  and  the  diurnal  variations  in  the 
current  systems.  In  both  figures  four  calculation  times  are  illustrated: 
thus,  for  example,  in  Figure  3-5a,  the  calculation  time  is  UT  (Universal  Time)  = 
0.0  hours,  summer  solstice,  when  the  Greenwich  meridian  (positive  x-direction) 
is  at  midnight.  In  Figure  3-5b,  the  calculation  is  at  UT  =  6.0  hours  when 
the  Greenwich  meridian  is  at  dawn  (i.e.,  the  sun  is  in  the  positive  y-di rection) . 
Similarly,  Figures  3-5c  and  3-5d  are  at  UT  =  12.0  and  18.0  hours  when  the 
Greenwich  meridian  is  at  noon  (the  sun  is  in  the  positive  x-direction)  and 
dusk  (the  sun  is  in  the  negative  y-direction) ,  respectively. 

In  Figures  3-5a  through  3-5d,  a  uniform  (constant)  potential  surface 
at  R  =  1  R£  was  used  as  the  boundary  condition;  while  in  Figures  3-6a  through 
3-6d  rotating,  uniformly  magnetized  sphere  was  used  at  the  boundary  of  R  =  1  R^. 
The  sensitivity  to  boundary  conditions  is  clearly  evident.  The  asymmetry  in 
both  sets  of  figures  is  due  to  the  induced  electric  field  resulting  from  time 
variation  of  the  magnetospheric  current  system.  This  is  particularly  notice¬ 
able  in  Figures  3-6a  through  3-6b  where  without  the  current  system  the  contours 
would  be  circular. 

In  Figures  3-5a  through  3-5d  the  change  in  potential,  <$> ,  between 
the  ionosphere  (where  it  is  assumed  to  be  zero)  and  the  equatorial  plane  is 
illustrated.  Comparison  of  these  Figures  3-5  and  3-6  suggests  that  if  the 
ionospheric  electric  field  (and  4>)  is  variable,  the  total  electric  field  in 
the  magnetosphere  will  also  vary  substantially  with  time.  Thus,  through  this 
mechanism,  the  ionosphere  may  exert  some  influence  on  the  dynamics  of 
magnetospheric  fields  and  low  energy  plasma. 
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Also,  Figure  3-5  indicates  the  change  in  <\>  between  the  ionosphere 
and  the  equatorial  plane  when  the  ionospheric  value  of  is  taken  to  be 
zero.  Thus,  magnetic  field  lines  are  not  equipotentials  even  when  only  the 
"small"  induced  fields  like  the  wobbling  dipole  are  present. 

The  various  components  of  the  total  electric  field  which  would  be 
"felt"  by  a  charged  particle  rotating  with  the  earth  in  geosynchronous  orbit 
are  illustrated  in  Figures  3-7  and  3-8.  The  uniformally  magnetized,  rotating 
sphere  boundary  condition  was  used  for  these  calculations.  Figure  3-7  dis¬ 
plays  the  field  as  a  function  of  longitude;  Figure  3-8,  the  electric  field 
as  a  function  of  time.  Contours  of  constant  I Etotai I  in  the  equatorial  plane 
are  shown  in  Figures  3-9  and  3-10  for  the  rotating  magnetized  sphere  boundary 
condition.  Figures  3-9a  through  3-9d  are  for  the  total  electric  field  in  an 
inertial  frame  at  four  times  (UT  =  0,  6,  12,  and  18  hours,  respectively).  The 
contours,  like  those  for  the  potentials  (Figures  3-5  and  3-6),  are  plotted  in 
geographic  coordinates;  all  comments  about  the  location  of  the  axis  with 
respect  to  the  sun  made  above,  are  applicable  here.  Figures  3-10a  through 
3-10d  are  for  the  total  field  in  the  rotating  coordinate  system. 

To  indicate  how  a  field  such  as  the  one  illustrated  in  Figures  3-9 
and  3-10  would  contribute  to  charged  particle  acceleration,  the  time  rate  of 
change  of  a  proton's  kinetic  energy  was  calculated  using  the  formula  from 
Northrop  (1963) 


d<w> 

dt 


where  <w>  is  the  kinetic  energy  averaged  over  a  gyration  period,  M  is  the 
magnetic  moment  of  the  proton,  q  its  charge,  and  V  the  velocity  of  its 
guiding  center.  The  results  of  this  calculation  are  shown  in  Figures  3-11 
and  3-12.  Figure  3-11  shows  the  energy  gained  by  the  guiding  center  of  motion 
as  a  function  of  the  guiding  center  kinetic  energy,  and  similarly,  Figure  3-12 
shows  the  energy  change  in  the  motion  about  the  guiding  center. 
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Figure  3-8.  Total  Electric  Field  Experienced  by  a  Charged  Particle  Rotating 
with  Earth  at  Synchronous  Orbit 
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Figure  3-9. 


I Etotal I  1n  a  rotat^n9  frame  for  a  rotating  uniformally  magnetized 

sphere  boundary  condition.  (Units  are  vol ts/meter  and  the  major 
divisions  are  5  R^.) 


33 


at i ng  magnetized 


3-11.  Time  Rate  of  Change  of  Proton  Kinetic  Enerny,  Guiding  Center  Motion  Contribution 
(UT  =  16.6  Hrs.,  Day  =  172,  Equatorial  Plane,  Longitude  =  0°). 


ENERGY  OF  ROTATION  ABOUT  GUIDING  CENTER  (eV) 

Time  Rate  of  Chanqe  of  Proton  Kinetic  Enerqy  Due  to  Rotation  About  the  Guiding  Center 
(UT  =  16.6,  Day  =  172,  Equatorial  Plane,  Ldnqitude  =  0°). 


3.3  EXTENSIONS  OF  THE  MODELS 

In  the  previous  section  the  electric  field  from  the  time  varying 
magnetic  field  resulting  from  the  diurnally  wobbling  dipole  was  discussed. 

The  procedures  and  methods  used  in  that  section  are  applicable  to  any  electric 
field  produced  by  a  time  varying  magnetic  field. 

Recently  (Olson,  et  al . ,  1979)  a  magnetic  field  model  was  developed 
that  represents  the  earth's  external  magnetic  field  during  disturbed  times. 

In  the  new  model  the  strengths  of  the  three  current  systems  can  be  varied 
separately.  It  is  thus  possible  to  descr  ioe  the  time  varying  magnetic  field 
associated  with  dynamic  processes  in  the  magnetosphere  such  as  magnetic  storms 
and  magnetospheric  substorms.  The  magnetic  vector  potential  is  also  deter¬ 
mined  and  can  be  used  to  describe  the  induced  electric  field  associated  with 
time  variations  in  the  magnetic  field.  The  electric  fields  associated  with 
the  storm  time  variations  are  considerably  larger  than  the  fields  described 
in  the  previous  sections. 

For  high  latitudes  the  field  line  topology  changes  dramatically 
with  small  changes  in  the  magnetospheric  tail  current  systems.  Figure  3-13 
shows  a  field  line  originating  at  72°  magnetic  latitude  for  varying  strengths 
of  the  tail  field  (0.7  to  1.4  of  normal  in  steps  of  0.1).  The  equatorial 
crossing  point  varies  from  about  13  R^-  to  over  60  R^. 

The  procedures  developed  for  modeling  the  quiet  time  tilt  dependent 
and  the  disturbed  time  magnetic  and  induced  electric  fields  dre  directly 
applicable  to  the  DNA  efforts.  The  electric  field  from  these  time  varying 
currents  are  found  to  be  important  for  plume  motions  and  plume  structures  in 
certain  altitude  and  latitude  regimes.  Furthermore  the  computer  codes  developed 
to  model  the  fields  are  directly  usable  for  calculating  the  magnetic  and 
self  induction  fields  from  a  plume.  Such  a  sample  calculation  is  given  in 
Section  5  for  a  sample  4  megaton  burst. 
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Section  4 
INTERACTION  OF 

ARTIFICIALLY  PRODUCED  PLASMA  WITH  NATURAL  FIELDS 


4.1  STRUCTURE  AND  MOTION  OF  ARTIFICIAL  PLASMA  IN  THE  NATURAL  FIELD 

4.1.1  Plasma  Characteristics  in  a  Low  Frequency  Field 

Electric  fields  observed  in  the  magnetosphere  have  been  discussed 
in  Section  2  and  techniques  which  can  be  used  to  model  these  fields  were 
discussed  in  Section  3.2.  Interaction  of  these  fields  with  artificially 
produced  plasmas  may  influence  the  structure  and  motion  of  these  plasmas 
even  though  the  internal  fields  of  the  plasmas  themselves  may  be  larger  than 
those  external  fields.  The  induced  electric  field  which  results  from  the 
diurnal  variations  in  the  magnetosphere  has  been  examined  quantitatively  as 
an  example.  This  field,  although  small  in  magnitude  0.1  mVolt/m),  is 
always  present  and  because  of  its  very  low  frequency  may  penetrate  artifi¬ 
cially  injected  plasmas  in  the  ionosphere  and  magnetosphere. 

The  degree  of  penetration  and  internal  magnitude  of  these  fields 
depends  on  the  detailed  structure  of  the  plasma.  It  is  shown,  however,  that 
under  conditions,  which  an  artificially  produced  plasma  at  high  altitudes  is 
sure  to  attain,  these  natural  electric  fields  can  become  important.  The 
regime  in  question  (called  "electromagnetic"  regime)  is  that  state  of  the 
plasma  in  which  the  mean  free  path  of  the  plasma  particles  and  the  plasma 
scale  size  are  much  greater  than  the  particle  Larmor  radius. 

All  naturally  occurring  gross  magnetospheric  electric  fields  are 
slowly  varying  in  time  when  compared  to  such  quantities  as  the  cyclotron  and 
collision  frequencies  of  the  plasma  particles.  Because  of  this  slow  variation 
of  the  field,  when  considering  the  interaction  of  the  background  field  and 
injected  plasma,  the  plasma  may  be  characterized  by  a  dielectric  constant  and 
the  degree  of  penetration  of  the  field  by  the  skin  depth.  The  resulting 


39 


electromagnetic  field  equations,  plasma  conductivities  and  dielectric 
constants  are  discussed  in  Appendix  A. 

To  simplify  the  problem,  certain  convenient  approximations  were 
employed;  the  general  procedure  is  that  followed  in  Frances  and  Perkins  (1975) 
Langevin's  equations,  which  are  adequate  approximations  for  the 
cold  and/or  long  wavelength  plasma  regime  of  interest  here,  are  used  to 
represent  the  motions  of  the  electrons  (subscript  e)  and  the  positive  ions 
(subscript  i) , 


dVe  _  _  _  _ 

'eUT  +  mevenVe  +  Ve1<W  '  -e<E+V6> 


and 


dVi 

"rdT  +  M1vinVi  +  Miuie<W  -  e(E+ViXB) 


The  neutrals  are  assumed  to  be  at  rest  in  the  coordinate  system  used 
To  obtain  the  conductivities  from  these  equations,  one  uses  Ohm's  Law 

J  =  a*E 

and  J  =  N(eV’--e7  ) 

i  e 

where  N  =  Ng  =  the  number  densities  of  the  ions  and  electrons. 
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To  further  simplfy  these  equations,  it  is  assumed  that 


_d_  _  3 

dt  ~  at 


+ 


v*v 


_3_ 

3t 


and  that  the  time  dependence  is  e’1wt  (i.e.,  consider  the  electric  field  in 

terms  of  its  Fourier  components),  where  to  is  of  the  order  of  2tt/(3600x24)  ^ 
-5  -1 

7.3  x  10  sec  .  Thus,  the  real  component  of  the  Pedersen  conductivity  is 
given  by 
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And,  similarly,  the  Hall  components  of  the  conductivity  are  given  by 
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where  use  was  made  of  the  fact  that  w  is  much  less  than  the  other  quantities 
given  above. 

The  collision  frequencies  appropriate  for  use  in  these  equations  are 
vgn  =  collision  frequency  of  electrons  with  the  neutrals 

=  2.6(104)R2NnTe1/2 

and  R  =  classical  ion  radius  (n,  lo’^m). 
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For  the  collision  frequency  of  ions  with  neutrals,  we  use 


vin 


=  6-K102)o2 
/"W- 


NnTi 


1/2 


where  W  is  molecular  weight  of  ions. 


The  collision  frequency  of  electrons  due  to  ions 
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The  collision  frequency  of  ions  with  electrons 
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In  what  follows,  four  regimes  are  considered: 

-5 

1.  The  low  altitude  (B  =  0.5  x  10  Tesla)  low  density  plasma 
(N  =  1012  m"3), 

2.  The  low  altitude,  high  density  plasma  (N  =  1014  m~3), 

3.  The  high  altitude  (B  =  1.0  x  lO-6  Tesla)  low  density  plasma 
(again  N  =  1012  m"3),  and 

4.  The  high  altitude,  high  density  (N  =  1014  m'3)  plasma. 

Further,  the  electron  temperature  is  taken  to  be  1/2  eV  or  5500°K 

18  -3 

and  the  neutral  density  in  the  low  altitude  regime  is  set  at  Nn  %  10  m 
and  Nn  ^  0  in  the  high  altitude  regime.  The  dominant  species  at  both  alti¬ 
tudes  is  assumed  to  be  atomic  oxygen.  The  resulting  values  of  the  relevant 
quantities  in  the  four  regimes  are  given  in  Table  4-1. 


The  most  important  parameter  for  considerations  of  field  penetration 
into  the  plasma  is  the  "skin  depth"  which  is  mathematically  defined  in 
Appendix  A.  It  is  that  distance  over  which  an  oscillating  electric  field 
loses  a  factor  of  e  of  its  energy.  As  can  be  seen,  the  electric  field  is 
expected  to  penetrate  well  into  the  plasma  especially  in  the  high  altitude 
regime. 

The  strikingly  large  skin  depths  apply  only  to  electromagnetic 
waves  of  the  low  frequencies  (7.3  x  10"^  sec'^)  considered  here.  For  fre¬ 
quencies  of  the  order  of  the  cyclotron  frequencies  the  character  of  the  whole 
problem  changes.  For  substorms  which  have  typical  frequencies  a  factor  of 
a  hundred  larger  than  those  above,  the  skin  depths  are  only  a  factor  of  ten 
lower  than  those  in  the  table.  These  skin  depths  are  also  large  and  similar 
conclusions  about  plasma  penetration  and  structure  can  be  made  about  those 
fields. 

4.1.2  Plasma  Behavior  in  a  Low  Frequency  Field 

From  simple  arguments  concerning  conditions  at  the  boundary  of  a 
plasma,  it  is  clear  that  the  component  of  the  total  electric  field,  ET, 
tangent  to  the  boundary  is  continuous  across  the  boundary  and  that  the  com¬ 
ponent  of  the  electric  displacement,  D,  perpendicular  to  the  boundary  is  also 
continuous  across  it.  From  such  arguments  it  can  be  shown  that  a  portion  of 
the  total  field  penetrates  the  plasma  in  an  anisotropic  manner  causing  dis¬ 
tortion  of  the  region.  In  the  "equatorial  region"  of  the  plasma  this 
penetration  becomes  important  because  the  resulting  ExF  drift  are  significant 
(relative  to  the  magnitude  and  direction  of  the  pressure  gradient  drift). 
(Equatorial  regions  are  those  regions  of  the  farthest  extent  of  a  given  mag¬ 
netic  field  line  in  the  vicinity  of  the  magnetic  equator.)  The  ExEf  drift 
produces  translation  of  the  plasma  and  may  also  cause  it  to  structure.  At  low 
altitudes  the  consequences  of  the  anisotropic  penetration  of  the  electric 

field  into  the  plasma  is  probably  not  significant  because  of  the  small  value 
_  o 

of  the  ExB/B  drift  velocities  there.  However,  the  significant  drift 
velocities  that  occur  for  the  plasma  in  the  equatorial  region  can  appreciably 
translate  the  locus  of  the  plume  during  late  times. 


Regions  of  sharp  gradients  (of  plasma  density)  of  a  finite, 
inhomogeneous  plasma  can  be  approximated  by  a  surface.  Using  Equation  A-2 
(see  Appendix  A)  and  Stoke' s  theorem  gives 

(f  E’d  =  /  (VxF)  •  nda 
s 

=  -ico  if  B*nda 


Since  B*n  is  continuous  across  the  surface,  so  is  the  tangential 
component  of  F.  Using  Equation  A-3  and  Gauss's  law 

(JsD-nda  =  J  V*Ddv 
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If  the  plasma  is  sufficiently  dense,  the  normal  component  of  F 
inside  the  plasma  is  much  less  than  the  normal  component  outside  the  plasma. 
The  tangential  component,  however,  does  penetrate  the  plasma  (see  Appendix  A 
for  skin  depth  estimates)  and  for  this  component,  geometrical  considerations 
become  important  (see  Figure  4-1).  Due  to  geometric  considerations  alone, 
the  tangential  component  of  F  will  be  a  maximum  in  regions  where  the  surface 
normal  becomes  perpendicular  to  F  and  a  minimum  when  it  becomes  parallel  to 
F.  If  the  skin  depth  (^  l//wy0/c)  is  large,  there  will  only  be  slight 
attenuation  of  this  component  in  passing  through  the  plasma. 
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In  those  regions  where  the  tangential  component  of  the  E  field 
penetrations  the  plasma  appreciable  distances,  the  drift  velocity  of  the 
plasma  due  to  this  field  may  become  significant.  The  equations  of  motion  of 
the  plasma  become 
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where  all  symbols  are  defined  in  Appendix  B. 


Figure  4-1.  Geometrical  Dependence  of  Penetrating  Tangential 
Component  of  E 
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In  Appendix  B  these  equations  are  solved  by  making  the  following 
assumptions: 


linEU  o 

at 


,  reasonable  for  slowly  varying  fields. 


v  •  ,  v 
in’  en 


,  are  negligible  in  the  equatorial  region  of  the  plume 


V  ~  v.  ,  reasonable  for  small  charge  displacement, 

“e  i 


ne  “  ni 


,  quasi  neutrality 


and  Te  »  Ti 


,  near  ambient  conditions. 


Under  these  conditions  the  solution  for  the  component  of  V 
perpendicular  to  B  is 


kT  vnxB 

V I  =  - 5 -  +  — ijr 

qB^ 


where  q  is  the  charge  of  the  particle  (+  e). 
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The  second  term  on  the  right  (ExB  drift)  is  independent  of  the 
particle  charge  and  is  the  same  for  all  species.  On  the  other  hand,  the 
first  term  on  the  right  (the  density  gradient  motion)  does  depend  on  particle 
charge  and  is  in  opposite  directions  for  ions  and  electrons.  Also,  since  Vn 
is  radial,  VnxB  is  azimuthal  and  does  not  contribute  to  translational  motions 
of  the  plasma.  The  penetrating  component  of  E  is  also  azimuthal  (see  Figure  4-1) 
hence  ExB/B2  is  radial.  The  two  types  of  drifts  (ExB  and  gradient)  are  not 
in  the  same  direction.  Thus,  even  for  "small"  E  fields,  the  ExET  drift  causes 
significant  displacement  of  the  plasma  over  long  times.  This  drift  is 
summarized  below  for  B  and  E  values  of  interest. 

ESTIMATES  OF  ExB/B2  DRIFT  VELOCITY 


Observed  Quiet  Time  E 

Range:  0.1-1  mV/m 

Value  Used:  0.1  mV/m 

Observed  Magnetic  Storm 
and  Substorm*  E 

Range:  1-5  mV/M 

Value  Used:  2  mV/m 

Ionosphere 

B  -  0.5xl0-4  Tesla 

2  m/sec 

40  m/sec 

High  Altitude  (L-3) 

B  =  1.1 xlO"6  Tesla 

• 

90  m/sec 

3 

1.8x10  m/sec 

*Enhanced  magnetic  activity  is  estimated  to  occur  about  25  percent  of  all 
time  near  solar  maximum  and  about  10  percent  of  the  time  nea *  solar 
minimum. 


How  significant  drifts  of  this  magnitude  are  depends  on  how  much 
time  the  plasma  spends  in  regions  where  the  ExB/B2  drifts  are  important. 

From  discussions  with  MRC  and  the  NRL  (National  Research  Laboratory)  groups, 
it  appears  that  a  specific  segment  of  the  plasma  will  spend  periods  of  time 
greater  than  an  hour  in  the  equatorial  region.  If  it  spends  just  one  hour. 
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it  will  drift  through  a  distance  1.8  km/sec  x  3.6  x  10  sec  =  6480  km  = 

1.02  R^.  This  corresponds  to  an  angular  translation  of  the  order  of  1.02/3 
radians  =  19.5°.  By  comparison,  neutral  winds  in  the  ionosphere  generate 
angular  drifts  of  the  order  of  l°/hour.  Whether  the  translation  is  in  lati¬ 
tude  or  longitude  is  determined  by  the  direction  of  the  electric  field. 

In  addition  to  the  translational  motions  of  the  plasma,  the  FxF 
drift  can  contribute  to  plasma  distortions  which  may  contribute  to  subsequent 
plasma  structuring.  This  is  illustrated  in  Figure  4-2  where  the  variation  of 
the  penetrating  component  of  F  results  in  the  varying  drift  pattern. 

In  summary  the  component  of  electric  field  tangent  to  the  plasma 
boundary  penetrates  the  plasma.  The  depth  of  penetration  is  most  significant 
at  higher  altitudes  and  varies  with  position  along  the  boundary.  In  the 

electromagnetic  regime,  where  both  plasma  size  and  particle  mean  free  path 

- 2 

is  much  greater  than  the  Larmor  radius,  the  ExB/B  drift  is  important. 

Because  of  the  geometrical  dependence  of  F  (and  the  resulting  geometrical 
dependence  of  the  drift  velocity),  the  motion  resulting  from  this  drift  not 
only  translates  the  plasma  but  also  distorts  its  shape  (contributes  to 
internal  structuring). 

It  is  suggested  that  further  investigations  be  made  to  see  in 
detail  if  and  how  the  time  dependence  of  the  distortion  leads  to  the  onset 
of  structuring.  Further,  the  FxF  drift  will  be  larger  still  in  the  environment 
produced  by  multiple  nuclear  bursts.  The  implications  of  this  for  real  battle 
situations  on  the  structuring  of  the  plume  and  the  distortion  of  electro¬ 
magnetic  communications  should  be  explored. 

4.2  PERTURBATION  OF  THE  NATURAL  FIELD  BY  THE  ARTIFICIAL  PLASMA 

4.2.1  Magnetic  Signatures 

Ground  station  magnetograms  of  the  world  wide  magnetometer  network 
have  observed  magnetic  field  perturbations  persisting  to  late  times  following 
the  various  nuclear  detonations.  Magnetograms  have  been  obtained  for  the 
post  STARFISH  period  from  the  World  Data  Center  A  in  Boulder,  Colorado. 

Figure  4-3  shows  the  variation  of  the  magnetic  field  at  Honolulu  and  Guam 
following  the  STARFISH  explosion.  The  zero  level  represents  the  level  of  the 
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(a)  Inhomoqeneous  Velocity  Distribution  Resulting  from  Externally 
Produced  ExB  Drifts 


(b)  Velocity  Flow  Lines  Resultino  from  Externally  Produced 
ExB  Drifts 


Figure  4-2.  Geometrical  Dependence  of  IxF/B^ 
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Figure  4-3.  Observed  Starfish  Magnetic  Field. 


ambient  field  before  and  after  the  explosion.  One  notes  that  for  late  times 
(t  >  500  sec)  a  detectable  signal  persists.  The  magnetic  field  at  late 
times  is  caused  by  the  currents  flowing  within  the  plume.  Computer  codes 
which  model  the  morphology  of  the  plume  can  calculate  the  currents  flowing 
within  the  plume  and  these  currents  can  then  be  used  to  calculate  the  expected 
magnetic  field.  The  calculated  magnetic  field  for  the  model  STARFISH  deto¬ 
nation  should  then  match  the  data  obtained  following  the  STARFISH  burst. 

Thus,  the  magnetic  field  signature  offers  a  test  of  the  validity  of  the  burst 
model.  The  magnetic  field  data  can  therefore  provide  an  independent  check 
on  the  total  plume  model  although  it  cannot  diagnose  specific  problems  within 
the  plume  model. 

To  test  the  feasibility  of  calculating  the  magnetic  field  signature, 
MRC  was  asked  to  supply  MDAC  with  information  on  the  currents  flowing  within 
a  plume.  MRC  first  supplied  MDAC  with  plots  of  the  currents  flowing  with  a 
typical  4  megaton  plume,  Figures  4-4,  4-5,  and  4-6.  Analysis  of  the  plots 
found  them  to  be  considerably  more  complicated  than  expected.  However,  an 
approximation  using  wire  loops  at  constant  altitude  to  represent  the  currents, 
indicated  that  the  field  under  the  plume  as  calculated  from  the  currents 
would  be  on  the  order  of  10's  to  100's  of  gamma.  The  complexity  of  the 
current  patterns,  however,  indicated  that  the  wire  loops  provided  only  a  crude 
approximation.  Thus,  a  more  accurate  procedure  using  direct  volume  integra¬ 
tion  over  the  current  system  was  developed.  MRC  agreed  to  supply  MDAC  with 
a  computer  tape  containing  the  current  vectors  at  a  grid  of  points.  The 
method  for  calculating  the  magnetic  field  from  the  currents  listed  on  the  MRC 
computer  tape  are  given  below.  (The  format  of  the  MRC  computer  tape  is 
described  in  Appendix  C.) 

Given  a  current  density,  T,  the  magnetic  field  ET,  at  a  distance  F, 
from  the  observer  is 


[  Jxr 


r 


dV 


(Note  all  equations  are  in  MKS  units.) 


(4-1) 
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E=-3  OE-IZ 

1=  1  OE-13 

M=  1-OE-ll 

10 

Fz-1  OE-12 

Jz  3  OE-13 

Nz  3  OE-11 

11 

0= -3  OE-13 

K=  l- OE-12 

0=  1  OE-IG 

11 

Hz-1  OE-13 

Lz  3.0E-12 

Pz  3.0E-10 

Figure  4-4.  Currents  Flowing  with  a  Typical  4  Megaton  Plume 
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Figure  4-5.  Currents  Flowing  with  a  Typical  4  Megaton  Plume. 
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Figure  4-6.  Currents  Flowing  with  a  Typical  4  Megaton  Plume 


elements  AV.. 

For  each  volume  element,  AV^ ,  we  associate  a  current  density  J-, 
where  7.  is  the  current  density  at  the  center  of  the  volume  element  and 
7. j  is  the  distance  to  the  center  of  the  volume  element.  Thus, 


_  7  N  J.xr. 

B  =  10’7  l  -L>  AV. 

i=>  I ri ! 3  ’ 


(4-2) 


A  similar  expression  is  available  for  the  magnetic  vector  potential 


7  N  T. 

fl  =  ,0’  J,  4Vi 


(4-3) 


MRC  supplied  the  current  vectors  in  the  a,  B  and  <p  coordinate 
system  with  the  currents  spaced  uniformly  in  a,  B  and  <j>.  The  a  coordinate 
is  along  a  magnetic  field  line,  the  £  coordinate  perpendicular  to  the  magnetic 
field  and  in  the  plane  containing  the  field  line,  $  is  perpendicular  to  a  and 
3.  < p  is  the  standard  spherical  coordinate,  a  and  3  are  defined  in  terms  of 

the  spherical  coordinates  R  and  0,  where  R  is  the  radius  vector  and  0  is  the 
colatitude.  The  following  equations  relate  a  and  3  to  R  and  G 


cosa  =  —k  cose 


(4-4) 


sinB 


-p-  sine 


(4-5) 


57 


Rj:  is  the  radius  of  the  earth  (6371  km). 

The  quickest  and  simplest  method  of  evaluating  the  integral  for  the 
Biot  Savart  law  (Equations  4-2  and  4-3)  is  to  perform  the  integral  in 
cartesian  coordinates.  The  volume  element  in  cubic  meters  is  given  by  the 
expression 


AV  = 


a4>r: 


sin  6, 


-  COS0  +  COS^G  -  ^  COS^G  +  y  COS^G 


*  '  I  -se  ’  3 


2  fcosG  ) 


5  (cosc^ 


3/2  i  S4 
cose 

0O 


(4-6) 


+  _L -  [  C0SQ  _  cos'-'G  +  £  cos'e  -  T  COS 

sin6e2  l 


-3a  ^  |  cos5e  - 


’•  1, 


The  symbols  are  defined  in  Figure  4-7.  A  complete  derivation  is 

given  in  Appendix  D.  Once  given  the  volume  element,  the  Ja  and  components 

must  first  be  transformed  into  cartesian  coordinates.  This  is  most  easily 

accomplished  by  first  transforming  and  0^  into  spherical  coordinates  J0 

and  J„. 
r 

Ja  can  be  split  into  a  component  in  the  r  direction,  Jar,  and  a 
component  in  the  0  direction,  JQ.  Since  a  line  in  the  a  direction  has  the 
equation 


R  = 


d  sin  6 

Ke  -  2 
sin^Bj 


(4-7) 
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and  since 


(4-8) 


it  can  be  shown  that 


i  —  2  j 

ar  sin9  ua0 


(4-9) 


Furthermore,  since 


J  =  J  2  +  J  2 

a  a  r  a  0 


(4-10) 


Substituting  Equation  (4-9)  into  Equation  (4-10)  and  solving  for  J  6  one 
gets  that 


CX0  j - a 

/  1  +  3cos  0 


(4-11) 


2cos0 


1  +  3cos 


A  similar  procedure  for  the  8  component  J.  gives  that 

P 


(4-12) 


2cose 


1  +  3cos  0 


T  Vb 


(4-13) 


and 


J 


6r 


-sine 


/l  + 


3cos  0 


(4-14) 


Since 


J  =  J  + 
r  ar  gr 


and  Je  =  Jc.e  +  Jee 


(4-15) 


One  can  combine  Equations  (4-12)  and  (4-14)  and  combining  Equations 
(4-11)  and  (4-13)  to  give 


Jr  * 


2  cos  0  Ja-  sin  6  Jg 


'rT~+ 


1  +  3  cos  8 


(4-16) 


J0  = 


2  cos  6  J&  +  sin  0  JQ 


•^77 


1  +  3  cos  0 


(4-17) 


The  final  transformation  to  cartesian  coordinates  is 


J  =  cos<{i  (sine  J  +  cose  Jj  -  sin<j>  J, 
x  r  0  cp 


=  si n<j>  (sine  Jr  +  cose  J^)  +  cos<> 


J,  =  cose  J  -  sine  Ja 
z  r  v 


(4-18) 

(4-19) 

(4-20) 


The  components  of  currents,  given  by  MRC  on  the  magnetic 

tapes  are  not  all  given  at  the  center  of  the  volume  element.  Instead  some 


61 


are  given  at  the  center  of  the  faces  of  the  cell.  Appropriate  interpolation 
schemes  are  used  to  calculate  and  at  the  center  of  the  volume 

element.  These  volume  centered  coordinates  are  then  converted  to  cartesian 
coordinates  using  Equations  (4-16)  through  (4-20).  These  currents  along  with 
the  element  calculated  using  Equation  (4-6)  are  used  in  Equations  (4-2)  and 
(4-3)  to  calculate  A  and  B  at  a  given  point  and  time. 

Table  4-2  lists  some  of  the  results  from  the  current  integration. 

All  values  are  along  the  4>  =  0  plane.  The  magnitudes  of  the  calculated 
fields  are  very  similar  to  the  STARFISH  data  set.  An  exact  comparison  cannot 
be  made  since  the  calculated  fields  are  for  the  model  4  megaton  burst  and 
the  data  are  for  STARFISH.  Furthermore,  the  effect  of  the  induced  telluric 
currents  have  not  been  included.  Also  note  that  the  burst  data  were  available 
only  at  500  sec  intervals.  It  is  therefore  not  possible  to  determine  if  the 
large  fluctuations  in  the  calculated  magnetic  field  values  (intuitively, 
monotonic  decay  might  be  expected)  are  "real"  or  if  there  is  a  problem  with 
the  codes. 


Table  4-2 

CALCULATED  MAGNETIC  FIELD  IN  GAMMA 
Typical  4  Megaton  Burst 


Colatitude 

Time 

Altitude 

500 

1000 

1500 

2000 

30 

0 

9.9 

5.7 

6.8 

13 

35 

0 

44 

27 

11 

36 

40 

0 

49 

102 

5 

60 

45 

0 

24 

42 

18 

34 

50 

0 

13 

19 

4 

16 

5 

0 

0.56 

0.52 

0.23 

0.89 

90 

0 

0.75 

0.35 

1.12 

2.47 

5 

30,000 

0.007 

0.012 

0.027 

0.035 

Maximum  Ej  =  9A/9t  =  0.2  mV/m 
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The  maximum  variation  of  the  vector  magnetic  potential  gives  a 
self  induction  electric  field  within  the  cloud  of  0.2  mV/m  which  is  on  the 
borderline  of  being  negligible. 

Thus,  a  method  has  been  developed  and  tested  which  permits  the 
calculations  of  the  magnetic  field  from  a  model  plume  at  any  magnetometer 
location.  These  calculations  can  be  used  to  check  the  magnetic  field  calcu¬ 
lated  by  the  computer  model  against  observational  data.  The  observational 
data  for  STARFISH  is  well  defined  and  thus  a  calculation  with  the  STARFISH 
currents  is  indicated  when  the  STARFISH  currents  become  available. 

4.2.2  Multiburst  Effects 

It  appears  clear  from  the  preceding  sections  that  interaction  of 
multiple  high  altitude  bursts  can  occur  through  mechanisms  other  than 
magnetohydrodynamic  waves  and  shocks.  High  altitude  fireballs  which  are 
largely  contained  by  the  magnetic  field  will  be  sensitive  to  any  changes  in 
the  magnetic  field  occurring  near  the  burst  point.  Remotely  occurring,  high 
altitude  bursts,  may  be  able  to  produce  sufficient  changes  in  the  magnetosphere 
to  directly  produce,  or  indirectly  precipitate  magnetic  field  fluctuations 
similar  to  the  naturally  occurring  magnetic  substorms. 

Large  displacements  of  the  magnetic  field  can  occur  in  the 
immediate  vicinity  of  a  high  altitude  burst  point.  Some  idea  of  the  magnitude 
of  this  displacement  can  be  obtained  by  comparing  the  energy  content  of  the 
ambient  field  with  energy  contant  of  the  ionized  atmosphere  and  bomb  debris. 
Depending  on  the  burst  altitude  and  weapon  yield,  displacement  of  the  order  of 
several  earth  radii  should  be  expected.  This,  coupled  with  the  magneto¬ 
hydrodynamic  wave  and  shock  effects,  can  yield  complicated  structure  and 
motion  in  existing  and  subsequent  fireballs. 

Any  changes  in  the  magnetosphere  influencing  the  structure  of 
natural  plasma,  in  the  current  systems  and  in  the  ionosphere,  will  have  an 
effect  on  both  earth  based  and  satellite  communication  systems.  To  properly 
evaluate  performance  of  such  systems,  the  importance  of  the  influence  of 
large  yield  multiple  bursts  in  the  magnetosphere  must  be  understood. 


Section  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  objectives  of  this  study  were  to  approach  the  problem  of  plasmas 
artificially  injected  into  the  high  altitude  ionosphere  from  the  viewpoint 
of  their  interaction  with  the  ambient  electromagnetic  field  environment 
(e.g.,  those  associated  with  the  magnetospheric  substorm  and  the  daily  wobble 
of  the  earth's  dipole  magnetic  field)  and  in  a  selfconsi stent  manner  to 
examine  the  mutual  influence  of  the  plasma  and  the  electromagnetic  field.  It 
was  found  that  the  natural  electric  fields  contribute  to  both  the  drift  motion 
and  structuring  of  the  plasma.  To  do  this,  the  principal  sources  of  the 
natural  electric  field  were  identified,  and  a  model  of  the  most  persistent 
field  (resulting  from  the  daily  variation  of  the  magnetospheric  magnetic  field) 
was  developed.  The  model  describes  the  total  electric  field  which  includes  a 
component  produced  by  the  plasma  as  it  responds  to  the  field  induced  by  the 
wobbling  dipole. 

The  approach  taken  in  deriving  this  model  is  outlined  in  the  proposal 
for  this  work.  The  model  was  made  sufficiently  flexible  to  allow  for  its 
expansion  to  include  contributions  from  other  sources  to  the  extent  that  they 
influence  the  structure  and  motion  of  the  various  magnetospheric  current 
systems.  Some  sources  of  the  natural  electric  fields  are  identified  in 
Table  2-1;  those  which  can  readily  be  modeled  with  our  approach  are  indicated 
in  the  table. 

The  natural  electric  fields  in  the  magnetosphere  have  periods  of 
the  order  of  0.25  hour,  or  more,  and  are  thus  very  low  frequency  in  character. 
As  a  result,  they  penetrate  the  plasmas  of  interest  here  with  very  little 
loss  of  energy  (i.e.,  have  large  skin  depths).  The  estimates  of  plasma 
properties  were  performed  under  several  simplifying  assumptions;  the  range  of 
values  used  are  typical  of  the  cases  of  interest  to  this  study.  The  dominant 
feature  of  the  problem  is  the  extremely  low  frequency  range  of  the  electro¬ 
magnetic  waves.  The  normal  variations  of  the  other  parameters  (such  as  the 
collision  frequencies)  with  respect  to  the  values  used  in  the  present  work  do 
not  alter  the  conclusions  appreciably. 
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In  addition  to  developing  the  electric  field  model  itself,  the 
possible  influence  of  this  electric  field  on  the  motion  and  structure  of  an 
artificially  generated  plasma  was  investigated.  Although  preliminary  in 
nature  the  results  of  these  calculations  indicate  that  the  plasma  can  be 
influenced  by  this  field;  the  ExB  drift  not  only  contributes  to  the  general 
motion  of  the  plasma  but  can  also  influence  internal  structuring  of  the  plasma. 
High  density  plasmas  in  the  ionosphere  are  expected  to  be  influenced  only 
slightly  by  the  smaller  natural  electric  field  there.  As  these  plasmas 
expand  outward  along  magnetic  field  lines,  the  penetration  of  the  electric 
field  increases  and  the  decrease  in  B  (i.e.,  ExB/B  increases)  produces  a 
significant  plasma  drift  velocity.  The  penetration  of  F  is  anisotropic  in 
nature  and  depends  on  the  detailed  structure  of  the  plasma.  The  plasma  is 
expected  to  drift  with  respect  to  the  magnetic  field  lines,  the  variation  in 
drift  velocities  throughout  this  region  contributing  to  the  internal  struc¬ 
turing  of  the  plasma. 

The  conditions  in  the  plasma  which  allow  penetration  of  the  external 
electric  field  are  also  those  under  which  the  assumption  that  magnetic  field 
lines  are  essentially  equipotentials  becomes  poor.  The  breakdown  of  this 
condition  allows  the  less  dense  plasmas,  in  the  extended  region  of  the  mag¬ 
netically  contained  plasma  plume,  to  act  more  or  less  independently  of  the 
denser,  low  altitude  regions.  Models  of  plasma  dynamics  in  the  magnetosphere 
which  depend  on  the  "field  line  coupling"  assumption  cannot  correctly  des¬ 
cribe  plasma  motions  in  the  regimes  of  interest  here. 

Although  not  included  in  the  original  work  statement,  work  was  also 
performed  on  the  perturbations  to  the  ambient  fields  resulting  from  the 
presence  of  the  artificial  plasma.  The  result  is  the  evaluation  of  the 
magnetic  signatures  which  could  be  observed  by  earth  based  magnetometer 
stations.  The  results  for  a  model  high  altitude  nuclear  weapon  burst  are 
compared  with  the  available  data.  These  calculations  are  based  on  current 
density  data  within  the  late  time  fireballs,  as  calculated  by  the  MRC  MICE 
and  MELT  codes.  The  results  are  compared  to  magnetometer  observations  during 
actual  bursts.  The  results  are  positive  and  indicate  that  this  technique  can 
be  useful  in  the  evaluation  of  the  performance  of  hydrodynamic  codes  such  as 
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MICE  and  MELT.  Further  extension  of  tnese  techniques  to  evaluate  the 
effects  of  high  altitude  nuclear  weapon  detonations  on  the  magnetosphere  may 
prove  very  useful  in  studies  of  communication,  satellite,  and  ballistic 
missile  defense  system  performance  under  nuclear  attack.  Multiple  burst 
effects  would  be  one  area  where  this  technique  could  be  applied.  Also  the 
techniques  used  for  quantitatively  determining  induced  electric  fields  can 
be  extended  for  use  in  the  auroral  regions  and  to  higher  latitudes  generally 
and  can  be  applied  to  the  multiburst  environment. 
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Appendix  A 

DERIVATION  OF  PLASMA  CHARACTERISTICS  IN  AN  EXTERNAL  MAGNETIC  FIELD 


r 


Maxwell's  equations  in  MKS  rationalized  units  are: 


VxH  =  J  +  e 


3E 
o  3t 


VxE  =  - 


9B 

3t 


VB  =  0 


(A-l) 

(A-2) 

(A-3) 

( A-4 ) 


Divergence  of  Equation  A-l,  using  Equation  A-3  gives 

V-0  +  =  0  (A-5) 

the  continuity  of  charge  equation. 

For  a  plasma  in  a  magnetic  field  and  an  alternating  low  frequency  electric  field 

0  =  oF 

and  p  =  -e  V- (XF) . 

Using  e1uit  time  dependence  in  Equation  A-5  gives 


icoe  X  =  a. 

o 


Let  a  =  °r  +  io-j  ,  X  =  X  +  iX^ 


then  X  =  — 

r  coe 


°1 


and 


X.  =  — 

1  toe. 


In  terms  of  Xr,  ar  one  has 


VxH  =  oy-E  +  iuje0(l  +  Xr)  E 


VxE  =  -iwB 


’•'S'1 


V'B  =  0. 


We  are  interested  in  low  frequency  u  «  v  «  u  and  long  wavelength  electric 
fields. 


Using  Langevin's  equation  (neutrals  are  assumed  approximately  at  rest) 


dV 

me  St  +  mevenVe  +  mevei(Ve  "  V  =  e(E+VexB) 


(A-5) 


dV, 


Mi  -ar  +  MivinVi  +  MiVVi  -  V  =  e<E+VixB> 


(A-7) 
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where  the  subscripts  e,  i  refer  to  electron,  and  positive  ion  properties, 
repsectively, 


ven  =  c0^1s10n  frequency  of  electrons  with  the  neutrals 

=  2.6(1 04 )  R2  N  T 1/2 
n  e 

R  -  classical  radius  of  the  neutrals  ^  10’10  meter 
and  Te  =  electron  temperature  (°K)  =  5,500°K  (1/2  e  v). 

Taking  =  10 12/cc  =  1818/m3,  vgn  =  2(104)  sec"^ .  For  plasma  in 

equatorial  regions  (L=3),  the  density  of  the  neutral  species  is  negligible  and 

v  -  0. 

en 

For  the  collision  frequency  of  ions  with  neutrals,  we  use 


6  1(102)R2 

.  =  ~  ~ ,v~ - 1 —  N  T.  ,  where  W  is  molecular  weight  of  ions. 


in 


/  W 


n  i 


Thus,  v-jn  =  1.1(102)  sec"^  in  the  ionosphere  (for  Nn  =  10^8/m3) 


and  v.jn  =  0  in  the  equatorial  regions. 


The  collision  frequency  of  electron  due  to  ions 


v  .  *  3.6(1 0-6 )  N.  T  -3/2  loge 

ei  l  e 


1.23(107)T3/2| 


17? 
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As  Te  -  5,500°K,  Ni  =  Ne  =  106/cc  =  1012/m3 


2  -1 

vei  ~  ^00  )  sec  in  tfie  boundary  regions  of  the  plasma.  In  the 

14  3  4  -1 

center  of  plasma,  =  Ng  -  10  /m  and  v  .  =  10  sec  . 


For  the  collision  frequency  of  ions  with  the  electrons, 
me 

vie  M  vei. 

-4  -5 

Now,  me/mp  =  5.44(10  )  and  me/m0+  =  3.37(10  ),  where  m  ,  mQ+  is  the  mass 

of  the  proton  and  oxygen  ion,  respectively. 

So,  for  the  low  density  plasma,  e.g.,  boundary  regions  of  the  plasma, 

Ni  -  N  =  1012  m  3,  v  ..  =  1.4(1 02 )  sec”1,  we  have  v-je(N+)  =  7.62(10’2)sec  and 

vie(°  )~  4.72(10  )sec.  For  the  high  density  plasma,  e.g.,  center  of  the 

plume,  N-  =  N  =  1014  m'3  and  v  .  -  104  sec’1,  we  have  v-  (H4)  =  5.44  sec”1, 

•  6  61  16 

vie(0+)=  0.337  sec"1. 


j 

Solving  Equations  A-6  and  A-7,  using  w  «  u>c,  v’s  and  3T-H  +  V-V;  a,so 
bB  pB 

letting  wce  =  >  w  ^  ,  gives  the  Pederson  components  of  the  conductivity 


n  (v  +  v  . ) 

/  \  \  _  e  '  en  ev 

xx  r  ‘  yy^  '  \  („  +  v  ,)5  + 17  ] 

'  en  ev  v  ce' 


2 

n.e  v.  (v-  +  v •  ) ( v  +  v  • )  +  ■  w 

i  ie  j-  v  in  ne/v^en  ev  ci  ce 

me  [( v  +  v  . ) 2  +  (jj2  ]  r  ( v  ■  +  v  •  ) 2  +  a)2 .  1 

LV  en  ev  ceJLV  in  ie;  ci  J 


A-4 


WWW*  wjpwwwj  u|i  i.ijij  mmm 


"ie 


(v •  +  v-  ) 

'  in  ie' 


m.  ,  v2  2 

i  (v.  +  v-J  + 


to. 


in  ie'  ci 


n  e  v  • 
e  ei 


mi 


(v-  +  v-  )(v  +  V  •)  +  (0  •  U)__ 

'in  ie/v  en  ei_  ci  ce 

2  ,  2  -  J 


C(ven  +  vei)2  +  “ce^^in  +  vie]"  +  wci] 


The  Hall  components  of  the  conductivity  are  given  by 


^a*y\ 


n  e 

( V>P  =  — 


ill. 


ce 


rn 


(v  +  V  •  )*"  +  u)2 
'  en  ei'  ce 


n,e2  vje 


0) 


(\J.  +  )  -  (\)  +  \1  •  ) 

i  ' n  v-ia'  » von  vai  ' 


ill. 


ce'  in  ie'  vven  ei'  “ci 
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[ (v  •  +  v  )2  +  to2  ][(v-  +  v-  )2  + 

L'  ei  en'  ceJ  '  in  ie'  “ciJ 


nie 


0), 


ci 


m. 


(v  •  +  v  •  )2  +  en2  • 

'in  ie'  ci 


u 


■  (v  +  v  . )  -  w  (v.  +  v  •  ) 

•1  '  on  m  »  'in  in' 


+  nee  vei  j-  ~ci  'Ten  'ei'  ~ce  "in  ’'ie 


m. 


i 


[(vei  +  ven)2  +  “ce^in  +  vie)2  +  wci  ] 
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In  the  ionosphere  B  ~  0.5(10  )  Tesla 


to  ~  9 ( 1 05 )  sec-1  ,  ui  •  (for  0+)  =  3(102)  sec  1  ,  w 

ce  ci  ^ 1 


i  ( H+ ) 


=  4.9(103)  sec"1. 
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In  the  equatorial  region  (L=3)  B  ~  1.1(10"^)  Tesla 


co  »  2(105)  sec'1,  w  .  (for  H+)  «  1 . 1  (102)  sec'1,  u  .(0+)  =  6.74  sec'1. 

C6  Cl  C  I 


In  the  ionosphere  v..  <  vgn.  Assuming  0+  species  only,  wC;.  ~  v^,  coce  »  v 

12  -3 

For  the  low  density  regime  (N  ~  10  m  },  ~  v^,  and  so 

coce  »  v0n  >  vin  ~  coci  ~  vei  »  vie  >  CO.  For  the  high  density  regime 
(N  ~  1014  rrf3),  vei  ~  vfin,  vie  <  vin,  and  so  coce  »  ven  -  vgi  >  vin  ~  u>c- 
>  vie  >  ^ec*ersen  conductivity  becomes  approximately 


Nie  <v1n  *  uie> 

- 

m. [w  .  +  ( v  •  +  v . _ )  ] 

iL  ci  in  ie'  J 


0  Low  Density 

N -  v.  1 ,05(10-3)  mho/m 

— —  •  ■  9  —  -y-  =  High  Density 

mi  v.  +  u>  .  1.05(10*1)  mho/m 

in  ci  '  ' 


In  the  equatorial  region  (high  altitude)  in  the  boundary  of  the  plasma  (N.  = 

12  3  -3 

10  /m  ,  the  Hall  conductivity  is  (o  )  "  1.4(10  )  mho/m.  In  denser  reaions 

xy  r 

14  3 

(N.  =10  /m  )  of  plasma,  the  Hall  conductivity  is  (o  )  ~  0.15  mho/r. 

•  xy  r 


The  Hall  conductivity  is  smaller  than  the  Pedersen  conductivity  for  the 
regimes  of  interest  here.  Thus,  in  the  ionosphere  the  skin  depth  parallel  tc 
the  electric  field,  which  is  perpendicular  to  the  magnetic  field,  is 

3 

/ — ^ —  3.3(10  )  km  at  plasma  boundary 

6  j  _  ip 

^o°xx  3.3(10^)  kn  in  dense  regions 

-5  -1 

where  we  used  w  =  2tt/3600 •  24  sec  =  7.3(10  )  sec  ,  as  the  characteristic  time 
dependence  for  the  E  field  at  quiet  times. 


A-6 


Thus,  the  electric  field  enters  well  into  the  boundary  layer  of  the  plasma. 


In  the  equatorial  regions,  for  L=3,  where  there  are  no  neutrals  and  where 

0+  ions  are  assumed  dominant  (v  -  v.  ~  0). 

en  in 

Low  Density  Plasma  Regime  uice  »  v  .  >  »  v1-e»a) 

High  Density  Plasma  Regime  w.  »  v  .  »ov,.;  ~  v.  »w 

3  3  ce  ei  ci  le 


n.e 

a(Pedersen)  = 


2,2 
v-  +  ai  • 
1  le  ci 


— 

*  1(10  )  mho/m  =  the  boundary  regions 

of  the  plasma 

=0.72  mho/m  =  the  dense  regions 


C 

1,1(10  )  km  at  the  boundary 
The  skin  depth,  6  =  « 

1.2(10  )  km  at  the  denser  plasma. 


This  was  done  using  ui  *  7.3(10  )  (24  hour  period)  as  the  characteristics 

time  dependence  of  the  E  field  during  magnetic  storms.  Again  the  electric  field 
enters  well  into  the  boundary  layer  of  the  plasma. 
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Appendix  B 

DERIVATION  OF  PLASMA  DRIFT  VELOCITIES 


To  obtain  the  expressions  for  the  significant  particle  drift  velocities,  we 
consider  the  following  equations  of  motion: 


3<"eV 


at 


+  v 


KT  n 
e  e 


-  —  (E  +  V  xB)  n  =  v. (V  -  V.)  nQ  +  v  ( V_ )  n„ 
m  e  e  ei  e  v  e  en  e  e 

e 


3(niVi) 


at 


+  v 


+  if-  (E  +  V,xB)  n  «  v.  (Vi  -  VJ  n  +  v.  V.  n, 
i  e  le  i  e'  e  in  i  e 


We  can  neglect  3/3t  terms  because  of  small  temporal  variation,  and  neglect  v 


ven  ™  tlie  ecluator''a1  regions. 


Assume  small  charge  displacement  ne  -  Vg  =  V.,  and  KTg  =  KT.. 


Thus,  for  each  species  KTg  +  q(E  +  VxB)  =  0,  where  q  =  +  e. 


Taking  the  cross  product  of  B  and  the  above,  gives 


KTe  ~^n^+  q  ExB  +  q  <VxB)xB  =  °- 


But  (VxB)xB  =  (VjxB)xB  =  VAB*  +  (Vi-B)B, 


where  is  component  of  V  in  plane  perpendicular  to  B. 


KT  VnxB/n  F  R 

Now  Vi  *B  =  0,  and  V2  =  • 

qB^  B^ 

Note  that  the  ExB  term  is  independent  of  the  sign  of  q  (i.e.,  common  drift  for 
each  species)  while  VnxB/qB  produces  a  drift  in  different  direction  for  each 
species. 

Since  Vn  is  radial,  VnxB  is  azimuthal;  ET  is  also  azimuthal  (high  at  top  and 
bottom,  zero  in  middle),  hence  ExB/B  is  radial.  That  is,  the  ExB/B  drift  is 
not  in  the  same  direction  as  VnxB;  and  as  a  result,  a  small  field  E  can  lead 
to  a  slow,  steady  effect  and  cause  significant  translational  motions  of  the 
plasma. 


(Note  that  ExB/B2  =  1.8  km/sec  for  E  =  2(10_3)V/m  and  B  =  1.1 (10'6)  Tesla.) 
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Appendix  C 

CURRENT  DENSITY  MAGNETIC  TAPE  FORMAT 


The  data  tape  of  current  vectors  supplied  to  MRC  is  an  unlabeled,  1600  BPI, 
binary  tape  written  by  the  CDC  Cyber  Record  manager  using  S  record  type 
and  C  type  blocking.  The  format  for  the  tape  is  given  below. 


C-l 


r 

i 


Record  1  -  Header  Record 

NCYCLE,  TIME,  TION,  NREC,  LENGTH1 ,  MREC,  LENGTH2,  JMXA,  JMXB, 
JMXP 

Record  2  -  Grid  Record 

(AKB(I),  1=0,  JMXA),  (BKB(I),  1=0,  JMXB),  (PKB(I),  1=0,  JMXP), 
(DAKB(I),  1=1,  JMXA),  (DBKB(I),  1=1,  JMXB),  (DPKB(I),  1=1,  JMXP) 
Stick  Records  ("NREC"  Records  of  "LENGTH1 "  words  each) 

(BKV(I),  1-1,4), (PKV(I),  1=1,4),  (AKB(I),  1=0,  JMXA), 

(CJA(I),  1=0, JMXA) 

Vertex  Records  ("MREC"  Records  of  "LENGTH2"  words  each) 

BV,  PV,  (AKC(I),  1=1,  JMXA),  (CJB(I),  1=1,  JMXA), 

(CJP(I),  1=1,  JMXA) 


! 


DEFINITIONS 


NCYCLE  MELT  master  cycle  number 

TIME  Simulation  time  (sec) 

TION  Ion  fluid  time  for  consistency  checks  in  MELT  (sec) 

NREC  Number  of  stick  records  -  JMXB*JMXP 

LENGTRI--. .^Number  of  60-bit  words  in  each  stick  record  =  10  +  2*JMXA 
MREC  Number  of-vectex  records  =  ( JMXB+1 )*( JMXP+1 ) 

LENGTH2  Number  of  60-bit  words  in  each  vertex  record  =  2  +  3*JMXA 

JMXA  Number  of  grid  cells  in  a-direction 

JMXB  Number  of  grid  cells  in  B-direction 

JMXP  Number  of  grid  cells  in  ^-direction 

AKB(I)  Array  of  cell  boundary  coordinates  in  a-direction.  Cell  "I" 
is  bounded  by  AKB(I-l)  and  ADB(I),  (Radians) 

BKB(I)  Array  of  cell  boundary  coordinates  in  B-direction.  Cell  "I" 
is  bounded  by  BKB(I-l)  and  BKB(I),  (Radians) 

PKB(I)  Array  of  cell  boundary  coordinates  in  <t>-di  recti  on.  Cell  "I" 
is  bounded  by  PKB(I-l)  and  PKB(I),  (Radians) 

DAKB(I)  Array  of  cell  widths  in  a-direction  (Radians) 

DBKB(I)  Array  of  cell  widths  in  (3-direction  (Radians) 

DPKB(I)  Array  of  cell  widths  in  ^-direction  (Radians) 

BKV(I),  Arrays  of  vertex  B  and  <p  -  coordinates  for  the  four  corners 
PKV(I)  of  each  stick.  (Radians) 

2 

CJA(I)  Array  of  Ja  current  density  components  (ABAMPS/cm  )•  These  density 

components  are  defined  at  the  points  (AKB(K),l/4  Z  BKV(I), 

I 


1/4  Z  PKV ( I ) ) 


BV,  PV  B  and  ^-coordinates  of  vertex  (Radians) 

AKC(I)  Array  of  cell  centered  et-coordi nates  (Radians) 


CJB(I), 

CJP(I) 


AKC (I ) = 1 /2 (AKB ( I )+AKB ( I - 1 )) 

Arrays  of  3  and  4>  current  density  components.  Defined  at 
points  (AKC(I),  BV,  PV).  (ABAMPS/cm2) 


I 
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Appendix  D 

DERIVATION  OF  LAGRANGEAN  VOLUME  ELEMENT 


The  derivation  of  Lagrangean  volume  element  (derivction  supplied  to  MDAC  by 
Bill  White  of  MRC)  is  given  below  (see  also  Figure  D-l): 

Volume  of  affo  Lagrangean  Element 


V  =  I  j  hahBh(j)  d«ded4>  = 
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~dG  sinG 
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;v  cosa^ 


p  /cosO 
ev  cosa,. 


dR  R  +  A<J> 


de  sine 


dR  R 


R 
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tosQ 


ev  cosa 
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d  /cose 
ev  cosa^ 

sine 

d  si n29 

Ke  •  2 
e  sin^B2 


A<t>  Re3  r°2 

--IT- j  dG 
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D-l 
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Figure  D-l.  The  (a, 3)  coordinates  for  the  field-aligned  dipole 
coordinate  system  far  from  the  earth. 
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